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ABSTRACT 

We study the escape of cosmic-ray protons accelerated at a supernova remnant 
(SNR). We are interested in their propagation in interstellar medium (ISM) after 
they leave the shock neighborhood where they are accelerated, but when they are 
still near the SNR with their energy density higher than that in the average ISM. 
Using Monte-Carlo simulations, we found that the cosmic-rays with energies of 
< TeV excite Alfven waves around the SNR on a scale of the SNR itself if the ISM 
is highly ionized. Thus, even if the cosmic-rays can leave the shock, scattering by 
the waves prevents them from moving further away from the SNR. The cosmic- 
rays form a slowly expanding cosmic-ray bubble, and they spend a long time 
around the SNR. This means that the cosmic-rays cannot actually escape from 
the SNR until a fairly late stage of the SNR evolution. This is consistent with 
some results of Fermi and H.E.S.S. observations. 

Subject headings: ISM: clouds — cosmic rays — ISM: supernova remnants 



Introduction 



Most of the cosmic-rays in the Galaxy are believed to be accelerated at the shock of 
supe rnova remnants (SNRs). Compared to cosmic-ray electrons that can be directly observed 
(e.g. iKoyama et al.l Il995l ). observational confirmation of accelerated protons is not easy. 
However, 7-ray observations have suggested that the protons illuminate molecular clouds 
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Recent 7-ray observations of molecular clouds around an SNR have given us information 
not only on the particle acceleration in the SNR but also on the escape of the particles from 
the SNR. By compa ring a simple model with latest Fermi and H.E.S.S. 7-ray observations, 
Fujita et al.l ( I2009at ) indicated that there are TeV cosmic-rays around the hidden SNR in 
the open cluster Westerlund 2, and the old SNR W 28. They showed that the diffusion 
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time-scale of cosmic-rays around the SNRs is much longer than that in the general region in 
the Galaxy. The d iffusion coefficient is less than ~ 1 % of that in the general region (see also 
Torres et al.l 120081 ). This suggests an important fact. The typical diffusion coefficient in the 
gener al region is D ~ f 29 cm 2 s _1 for particles with an energy of ~ 1 TeV ( IBerezinskii et al. 
19901 ). If the coefficient is 1% of that typical value or D ~ 10 27 cm 2 s -1 , the time-scale on 
which the particles cross the scale of an SNR (R s ^ ~ 15 pc) is ~ R^ h / (6D) ~ 1 x 10 4 yr. This 
is comparable to the time in which a shock with a velocity of ~ 10 3 kms _1 crosses that scale. 
This means that particles that can leave the shock neighborhood would not easily escape 
from the periphery of the SNR. Here, the shock neighborhood means the region where some 
non-li near effects generate strong ma gnetic waves or cause strong amplification of magnetic 
fields (ILucek &: Belli I200CH : iBelll 120041 ) . In that region, the particle diffusion would follow the 
so-called Bohm diffusion and the particles are efficiently accelerated. 

In this letter, we theoretically investigate particle diffusion in the ISM outside this shock 
neighborhood. We show that cosmic-ray streaming generates Alfven waves around the SNR, 
which scatter the particles and make the particle diffusion in the ISM around the SNR 
significantly slow. We also show that the cosmic-rays actually can not easily escape from 
the p eriphery of the SNR. In contrast wi th previous studies (e.g. iPtuskin fc Zirakashvili 
20031 ; see also iLee. Kamae &: Ellison! 120081 ). we consider the time-evolution of the diffusion 
coefficient and the effect of cosmic-ray diffusion in the ISM before being affected by the 
streaming. We use Monte-Carlo simulations to study the diffusion, while at the same time 
we calculate the evolution of the SNR and the excitation of Alfven waves based on a simple 
model. We consider the diffusion in the ISM around an SNR that are highly ionized, because 
ion neutral friction does not significantly affect the growth of the Alfven waves. Observed 
molecular cl o uds th at are bright in the 7-ray band may be immersed in the ionized ISM. In 
Fujita et al.l ( 12009al ). it was not specified whether the observed slow diffusion of cosmic-rays 
is happening inside or outside the molecular clouds. In this study, we take the latter stance. 



2. Models 

Cosmic-rays are scattered by Alfven waves generated by the streaming of the cosmic- 
rays themselves. We treat the propagation of a cosmic-ray particle as a three-dimensional 
random walk. We assume that the diffusion coefficient depends on the position r, although 
the scattering is isotropic at each point. For a given diffusion coefficient D(t, r, E), where E 
is the energy of the particle, the mean free path of a particle is written as l m = 6D/v, where 
v is the velocity of the particle. The diffusion coefficient is given by 

„ 4 pvc , . 

Sir eBtp 



-3 - 



where p is the momentum of the particle, c is the speed of light, e is the positron charge, B is 
the magnetic field, and ip(t, r, E) is the energy density of Alfven waves per unit logarithmic 
bandwidth (which are resonan t with particles of energy E) relative to the ambient magnetic 
energy density Um HBellj ll9T8h . At the rest frame and the outside of a shock, 

dip 
~dt 



a 



(2) 



where F is the damping rate (jBelllll978l ). The growing term is given by 
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where vx is the Alfven velocity, and f(t,r,E) is the distribution function of cosmic-rays 
(jSkillinglll975l ; lBelllll978l ). For the damping, we consider the one owing to collisions between 
charged and neutral particles; 
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where is the neutral hydrogen density, T is the temperature of the ISM, and 

-0.4 / n \ 2 



PcohC = 8 



_nu_ 

cm _: 



-3/2 



T 



10 4 K 



B 



1/xG 



GeV 



(4) 



(5) 



(IKulsrud fc Cesarskylll971t lO'C Drury et al.lll996l ). For p > p coh , neutrals participate in the 
coherent oscillations of the ions in the Alfven wave, which decreases the damping rate. 



The evolution of an SNR follows the Sedov solution for t sc d < t < t ra d; 



Rsh(t) — Rrad(t/t 
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where t se d is the starting time of the Sedov phase, and i? ra d is the radius of the SNR when 
its radiative cooling becomes effective at t — t ra d- For particle acceleration in the shock 
neighborhood, we adopt a simple model. We assume that on ly particles with a maximum 
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particle diffusion outside the shock, we assume that the growth of Alfven waves stops inside 
the shock for simplicity {ip = 0). 

We assume that ip = at t = 0. When ip is close to zero, the diffusion coefficient in 
equation ([T|) diverges. Thus, we use a diffusion coefficient of D = mm(D a , D s ) in actual 
calculations, where D s is the standard value in the Galaxy at positions far away from SNRs: 
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( IGabici et al.l 120091 ). We do not consider the cooling of particles and the reacceleration of 
particles that are caught up by the SNR. 



3. Results 

First, we consider cosmic-ray diffusion in the ISM that is completely ionized (r = 0). 
We assume that R r!U \ = 15 pc and t ra d = 1 x 10 4 yr, which are close to those for Westerlund 2 



and W 28 JFuiita et d]l2009afl ). The strength of undisturbed magnetic field is B = 10 /i G. 
The mass density of interstellar medium (ISM) is p = 8.35 x 10~ 24 g cm -3 . The maximum 
energy i? max decreases from 100 TeV (t = t scd = 215 yr) to 1 GeV (t = t ra d), which leads 
to 5 = 3. For the calculations of if) (equation [2]), we set the inner boundary at r = 3.2 pc, 
where the highest energy cosmic-rays are injected (E = 100 TeV). We consider the gradient 
of / only in the radial direction. The outer boundary is set at r = 50 pc. The innermost 
mesh has a width of ~ 0.2 pc, while the width of the outermost mesh is ~ 6 pc. The 
time-step for the calculation of ip is taken as the time in which E m3uX (t) decreases by a factor 
of 10 ' 2 , which is a fa ctor of a few smaller than the duration in which particles with i? max (t) 
are actually injected (jCaprioli. Blasi &: Amatll2009l ). 



In Fig. [lj we present the distribution of particles with an energy of E = 1 TeV at 
t = 2.2 x 10 3 and 1 x 10 4 yr (or t ra d)- The shock radius at those times is = 8.1, 
and 15 pc, respectively. The particles with that energy are injected at t = 1 x 10 3 yr at 
i? s h = 6.0 pc. The number of the injected particles with that energy is 100,000. Immediately 
after the injection at t — 1 x 10 3 yr, the particles rapidly leave the shock because of the large 
diffusion coefficient (D = D s ), until Alfven waves grow enough to scatter the particles. 

Figure [1] indicates that many of the particles are still wandering around the SNR at 
t = 2.2 x 10 3 yr. Until this time, the streaming of the cosmic-rays has generated Alfven 
waves outside the SNR, and they are strong enough to scatter particles (but still ^ < 1). 
As a result, the diffusion coefficient has become much smaller than D s in equation ([7]). In 
fact, at r ~ 10 pc, D a /D.^ ~ 0.01 for E = 1 TeV, which is consistent with the Fermi 



and H.E.S.S . result s (jAharonian et al.l 120071 . l2008t lAbdo et al.l |2009b() as is discussed in 



Fujita et al.l (j2009al ). Even at t = t ra d, a significant fraction of the particles remain around 
the expanding SNR. Because of the relatively high value of ip and the low value of the 
diffusion coefficient, cosmic-rays cannot easily escape from the region peripheral to the SNR. 
At t = t ra( j, 71% of the injected particles with E = 1 TeV are around the SNR (r < 20 pc). 
The number density does not decrease smoothly toward larger r (Fig.[T]), because the diffusion 
becomes slow where the gradient of / is large (equation [3]), which tends to steepen the 
gradient further. 
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Figure |2] shows the spectrum of cosmic-rays at t — £ ra d- We assume that the total energy 
of cosmic-rays accelerated for t sed < t < t rad is 10 50 erg. Although we give the same number of 
particles for a given energy in the calculations, we give them weights so that the distribution 
function of all the particles that are accelerated at the shock is fo oc E~ 2 . Figure [2] indicates 
that / oc E~ 2 - 75 for E > 10 TeV regardless of positions. The in dex (—2.75) is the one when 



particles haye alrea dy prevailed beyond the region considered (lAharonian fc Atoyanl 11996 



Fujita et al.l l2009al ). For this energy range, Alfven waves are not much excited, because 
equation (jTJ indicates that the diffusion coefficient when ip m (D — D s ) is large, and 
because the particles rapidly escape from the SNR without making a gradient of / large 
enough to excite Alfven waves (equation [3]). On the other hand, the diffusion coefficient 
when if) w is smaller for E < 3 TeV (equation [7]), Alfven waves are generated through 
the gradient of / in this energy range. As a result, cosmic-rays in this energy range are well 
scattered by the Alfven waves, and spread only slowly as a cosmic ray bubble. Thus, the 
spectrum follows the original form of / oc E~ 2 , and the value of E 2 f is distinctively large at 
r < -Rsh (= 15 pc) compared to the one at r > R^. The position of the cutoff in the spectra 
(~ 3 TeV) depends on the diffusion coefficient of the unperturbed ISM (equation [7]). This 
indicates the need to study the escape of cosmic-rays globally, not locally around the shock. 
At r = 32-36 pc, the spectrum is harder than E~ 2 at E < 3 TeV, because most of the low- 
energy particles have not reached that region. It is to be noted that the wave energy does 
not exceed the background magnetic energy around the SNR in this energy range {ip < 1). 
This is different from the shock neighborhood mentioned in § [H where if) > 1 is expected. 
The retained cosmic-rays may be released at t > t ra( j when the SNR can no longer hold them. 

The spectra are different when the ISM is not completely ionized. In Fig. [3j we assume 
that the ISM temperature is T = 10 4 K, and the neutral hydrogen density is wh = 0.1 cm -3 
(equations [I] and [5j). We chose these parameters because the neutral damping is marginally 
effective, and we can see the effect clearly. Other parameters are the same as those studied 
above. In Fig. EJ the neutral damping is effective for E < 1 TeV at r > 16 pc. Around the 
SNR or r ~ 12-16 pc, the generation of Alfven waves overwhelms the damping at E ~ 1 TeV 
because of the large gradient of /. At r ~ 12-16 pc, particles that are recently injected are 
seen at E < 10 GeV. 



4. Summary and Discussion 

We have shown that cosmic-rays with energies of < TeV excite Alfven waves in the ISM 
around an SNR even after they leave the shock neighborhood where they are accelerated. 
The particles are scattered by the waves, which makes the diffusion of the particles slow. 



-6 - 



As a result, the particles remain around the expanding SNR for a long time. This means 
that even if the particles are accelerated in an early stage of the SNR evolution, they cannot 
actually escape from the SNR until a fairly late stage. 

The ISM must have been well ionized for the cosmic-rays to be held around the SNR, 
because the neutral damping of Alfven waves is not effective. For Westerlund 2 and W 28, 
Fujita et al.l (j2009al ) indicated that the supernova exploded in a hot cavity, which had been 
created through activities of the progenitor star and/or explosions of other supernovae. 
Therefore, the ISM might have been ionized. While the maximu m energy of protons around 
the SNR is ~ 2.7 TeV for W 28, it is ~ 47 TeV for Westerlund 2 jFuiita et al. l l2009ah . which 
is larger than that predicted in Fig. [2] (~ 3 TeV). However, the maximum energy increases 
to > 10 TeV if the magnetic fields in the ISM is B ~ 20 fiG and the ISM is fully ionized. 

Recent H.E.S.S and Fer mi observations showed that 7-rays are produce d in molecular 
clouds adjacent to SNRs (e.g. lAharonian et al.ll2004j . 120081 ; lAbdo et al.ll2009al ). In molecular 
clouds, Alfven waves may totally dissappear through neutral damping. If this is the case, 
they may be illuminated by cosmic-rays in the ionized ISM surrounding the clouds. The 7- 
rays may be generated through pp-interaction in the clouds. The cosmic-rays that enter the 
clouds may quickly pass the clouds unless some other disturbers such as turbulence scatter 
the cosmic-rays. In our study, the cosmic-rays are confined around an SNR on a scale of the 
SNR itself (Fig. [[]). Thus, molecular clouds that are bright in the 7-ray band may not be 
found in the region beyond this scale. 

If cosmic-rays continue to be held by an SNR, they may be affected by adiabatic cooling 
as the SNR expands. Moreover, some of the cosmic-rays engulfed by the SNR would be 
reaccelerated. These may affect the cosmic-ray spectrum in the Galaxy. 
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Fig. 1. — Number density of particles with an energy of 1 TeV at t — 2.2 x 10 3 (solid) 
and 1 x 10 4 yr (dotted). The radius of the SNR at those times is R s h = 8.1 and 15 pc, 
respectively. The supernova exploded at the origin. 
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Fig. 2. — Energy spectra of cosmic-rays at t — t ra d = 1 x 10 4 yr. The ISM is fully ionized. 
The solid line is for r =12-16 pc, the dotted line is for 16-20 pc, the dashed line is for 
24-28 pc, and the dot-dashed line is for 32-36 pc. 



-lo- 



co 
I 



10°h 



o 

> 

CD 



£io- 2 



10 



-4 




12-16 pc 
1 6-20 pc 
24-28 pc 
32-36 pc 



10' 



10° 
E(MeV) 



10 



8 



Fig. 3. — Same as Fig. [2] but the ISM is not completely ionized. 



